Current induced decoherence in the multichannel Kondo problem 
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The properties of a local spin S = 1/2 coupled to K independent wires is studied in the presence 
of bias voltages which drive the system out of thermal equilibrium. For K ^> 1, a perturbative 
renormalization group approach is employed to construct the voltage dependent scaling function for 
the conductance and the T-matrix. In contrast to the single-channel case, the Kondo resonance is 
split even by bias voltages small compared to the Kondo temperature Tk, V <C Tk- Besides the 
applied voltage V, the current induced decoherence rate rcl 7 controls the physical properties of 
the system. While the presence of V changes the structure of the renormalization group considerably, 
decoherence turns out to be be very effective in prohibiting the flow towards new nonequilibrium 
fixed points even in variants of the Kondo model where currents are partially suppressed. 
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The problem of a local spin S exchange coupled to 
K independent conduction electron reservoirs, inspite of 
its apparent simplicity shows a rich variety of phases. 
For example when K=2S, the local spin is completely 
screened and the system behaves as a Fermi-liquid [l[ . In 
contrast for K > 2S, the spin is over-screened, and the 
system exhibits non-Fermi liquid behavior characterized 
by a zero temperature entropy S = In g where g is a non- 
integer [2-5]. Recently the over-screened Kondo problem 
for K~2, S—l/2 was realized in a controlled experimental 
set-up [6|, opening up the possibility of probing these 
exotic systems in the far out of equilibrium regime. 

An important question in the study of any strongly 
correlated system is the possibility of realizing new fixed 
points and scaling behavior by driving the system out 
of equilibrium. In [7:] it was proposed that the single 
channel Kondo model in the large bias limit should flow 
to a new fixed point which is characterized by a change 
in the number of independent screening channels as co- 
herent scattering processes between leads are prohibited 
by the bias voltage. However subsequent studies [§-[l2j 
of the nonequilibrium single-channel Kondo model ruled 
out such a nonequilibrium regime due to current induced 
decoherence which qualitatively plays the role of an effec- 
tive temperature (or an infra-red cutoff), and can be as 
large as the voltage in the single-channel Kondo model. 

In this paper we study the over-screened Kondo prob- 
lem when a local spin 5=1/2 is coupled to K independent 
current carrying channels. In the limit K 3> 1, the in- 
termediate coupling non-Fermi liquid fixed point can be 
studied within renormalized Kcldysh perturbation theory 
in 1 /K . We show that the current induced decoherence in 
this model is 0(1/ K) and thus considerably suppressed. 
In spite of this we find that decoherence is highly effec- 
tive in cutting off the flow to any nonequilibrium fixed 
points. We present results for the voltage dependence of 
the conductance and the T-matrix, the latter being re- 
lated to the local density of states that can be probed in 



a tunneling experiment. 
The Hamiltonian is H 



Hq + Vx where 



H = £ (e fc -AC)4 



tjma Ck <ymu 
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represents K independent free electron reservoirs labeled 
by to. Each of the K reservoirs is split into a left (L) 
and a right (R) part which can be maintained at dif- 
ferent chemical potentials p™ ^ p 7 ^. The local spin is 
coupled to the spin density of each wire via the exchange 
interaction, 



v x = \s- 



Jaf3C k(Tma d aa ,C k , a ' mli (2) 
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Here, Jlr connects the L and R leads so that a net cur- 
rent can flow within each reservoir when p™ ^ p 7 ^ . How- 
ever there is no flow of current from channel m to channel 

to' 7^ TO. 

We evaluate two physical quantities, the current 



through the m-th lead /„ 

I m =—i%JLRi> ■ J2 
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which is given by 
o<j<j>Ck><jiRm - ti.c.) , and 



the T-matrix of the electrons in the m-th lead defined 



by G 
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where aa is the retarded propagator for 

the electrons in the m-th lead and the T-matrix 
is the expectation value of the composite op- 
erator T™« Q '(t, t')=-iT c {O a (t)Ol(t')) where 

O a (t)= 2 5Zaifci<7i ^""l ' (J cra 1 c ma 1 k 1 a 1 a{t) j b i 

denote Keldysh labeling and Tq denotes Keldysh 
time-ordering 13]. Formally, under assumptions that 
the interactions Vx was switched on at t=— oo, the 
expectation value of any operator O at time t=0 is 13 1 



(O(0))=Tr Ufe^-oo dt ' v ^^ Q foe " \ p 

where p is the initial density matrix at t——oo and 
Vx(t)=e tHot Vxe~ tHot . We assume p=pi ea ds <8> Ps where 
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ps is the density matrix for the free spin, and pi ea ds is 



such that (c 



,t 



kcrma C k<?ma 

Fermi function. We also assume that the leads have a 
uniform density of states v and a bandwidth 2D. 

Identical voltage drops across the K-wires: Let 

us suppose that all the K wires have the same voltage 
drop V—\ pl ~ M.R I applied across them. Then a per- 
turbative treatment to two-loop order, where only the 
diagrams which are O(K) at two loop are kept gives the 
following result for the conductance G m — dI m /dV, 
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In addition the imaginary part of the T-matrix of the 
m-th wire is found to be, 



-7wlm 



[j7 a,<7 (n) 



3tt2 
16 



D 



I ^ - Mai I 



' K tl c(cx.i J a.-\ a' In 



Q - /i Q1 



L"v;e(D- i/^-w I) 

7(5 



(4) 



According to Eq. ([3]) all logarithmic singularities in 
the conductance are cutoff by the voltage. This reflects 
the fact that in a model where the cutoff D is much 
smaller than the voltage difference V no current will flow 
by energy conservation. In contrast, the calculation of 
the T-matrix shows that even in this regime resonant 
spin-flip processes lead to logarithmic renormalizations. 
For the renormalization group (RG) analysis it will be 
important, that some Jlr processes do not contribute for 
D < V as described by the 9 function terms in Eq. (|4|). 

We will formulate the RG equations in terms of dimen- 
sionless couplings, g a b=vJab assuming symmetric cou- 
plings, gLL=9RR=9d but we will allow for g LR ^ g d . 
While for a simple Anderson model one always obtains 
gLR—gdi this is not valid for more complex models. 

When deriving two-loop equations in full generality for 
a non-equilibrium situation using, e.g., functional renor- 



malization group [141 ] . flow equations [l(| or real-time 
renormalization group approaches [l2[ , it is necessary to 
take into account the full energy dependence of interac 



tion vertices. Even at the one- loop level it is useful [15 1 
to keep track of how the coupling constants g depend on 
the energy of the incoming electron. To leading order in 
X/K, however, and for the quantities considered in this 
paper, we believe that it is sufficient to use a simpler 



version guided by results from perturbation theory, con- 
centrating on a few on-shell coupling constants. Already 
in perturbation theory to order J 3 the structure of loga- 
rithmic corrections depends on which physical quantity is 
considered. The simplest case is the conductance where 
all logarithmic corrections up to the order considered are 
cutoff by the voltage. Thus for the conductance it is suffi- 
cient to use the well-known equilibrium RG equations [H 



dgd 
d\nD 

dgLR 
d\nD 



{9d + 9lr) 
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g d (2gj + 2g 2 LR ) 
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supplemented by the condition that the RG flow is cutoff 
at D=V. 

However, when one requires instead that the T-matrix 
T aa (p a ), evaluated at the Fermi energy of the lead, re- 
mains invariant under RG, one obtains the RG equations 



dgd 
d\nD 

dgLR 
dlnD 



K 

9lr®v) ~ -^9d (2ffd 
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Here the factors 8y = Q(D — V) take into account that 
part of the logarithms and therefore the RG-flow is cutoff 
by the voltage when the cutoff becomes smaller than the 
voltage. Therefore the equilibrium RG flow is modified 
at the scale V in a way which can be simply understood 
as a suppression of all resonant scattering processes from 
one to the other chemical potential by gLR- 

Note that Eq. ([6]) is not valid in the limit D — > be- 
cause the RG flow is ultimately cutoff by the decoherence 
processes as had been pointed out previously [8l-fllj|. A 
careful analysis of the decoherence rates (including 3 loop 
diagrams) has been performed by Schoeller and Reining- 
haus in [l2j by investigating the time evolution of the 
density matrix of the impurity. Their analysis shows that 
within the precision of a two-loop calculation (i.e. within 
our model to leading order in 1 / K) , the decoherence rate 
is given by the simple formula used previously [H E| for 
one-loop calculations. Up to prefactors of order 1, irrele- 
vant for our discussion, the RG flow is ultimately cutoff 
at the scale T, 



T = nKg L R(V) 2 V 



(7) 



where gLR,{V) is the running coupling constant at the 
scale D=V. The factor K arises as each of the K chan- 
nels contributes to the dephasing. V is the Korringa spin- 
relaxation rate [§- 11 and is proportional to the current. 

Let us make the assumption gd=gLR. (we will relax this 
condition later). For the regime D > V it is convenient 
to define g—2gd which obeys the RG equation 

dg 



dlnD 



K 3 
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FIG. 1: (Color online) Main panel: The conductance 
in the scaling limit and for several different A=2/K. 

Ga=(t£) ^A 2 . Inset: Im[T(n)] for V/T K =0.1 and the 

same A as in the main panel. Ta=37t 2 A 2 /32. 



The above equation has the well known fixed point at 
g*=2/K 0, whereas the scaling dimension of the leading 
irrelevant operator is A=/3' (g*)=^ valid for K 3> 1 (the 
exact result [4] is A = 2/{K + 2)). Integrating Eq. © 
upto an energy-scale D, and defining the Kondo temper- 
/ \ K / 2 

ature as T K =D (p) e~ 1/ao , (D ,g being the ini- 
tial bandwidth and coupling constant respectively), one 
gets 



\g{D)-g 

= l.9o - 9* 



\g(DY 


KA/2 
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Eq. © is valid for arbitrary 
TV, i.e. for g <C g*, we obtain g(-D) 

g \n(2/KlnD/T K ) , n/ 1 ^ 3 



For L> > 
l , 



\tlD/Tk 1 

_ (2/KlnD/T K ) 0( 1 n 3 p Q T setting 

2 (in d/Tk) 2 + t -AinD/T A - J ■ rot u ^ J -ff: setting 
g(D)=g* on the r.h.s of Eq. © one gets 



- .9* 1 



(10) 



where £ = (1 — go/ 9*) e~ A ^ 9 . In the scaling limit 50 ~> 
0, and for large K, £ takes the universal value £ = 1/e. 

Integrating upto the energy scale D—V, the loga- 
rithms in the conductance are resummed giving G m ps 
^T1T9lr( v )- In particular the result for V < T K is 



G„ 



3tt£^ 
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TV 
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where we have dropped higher order terms ~ ('Tk) 
Note that the conductance near the fixed point is a 
quantity of 0(1/ K 2 ). For arbitrary V/Tk and in the 
scaling limit, the conductance is given by the univer- 
sal function G/G A = l/[l + W {e- 1 {V/T K ) A )] 2 where 



G A =(pf) TT A2 and W{x) is the Lambert W function. 
The conductance for several different A is plotted in 
Fig. [TJ At the fixed point g* LR =\/K. Thus the deco- 
herence rate (defined in Eq. ©) for V <C Tk is, 



ttV 



(12) 



and thus apparently small for K 1. However, as we 
show below, this decoherence rate plays an important 
role in cutting off the logarithms in the T-matrix. 

For V=T=0 the T matrix has a powerlaw cusp, 
ImT(w) w d - c 2 |w| A Eqn. © suggests that this 

peak will split by V and we also expect it to be broadened 
by T. As r <C V this splitting will be observable even 
for V <C Tk- In the following we will use RG to calcu- 
late this split Kondo resonance for V <C Tk (results for 
V 3> Tk follow from the literature of the single-channel 
Kondo model [HI]). We assume the symmetrical applica- 
tion of voltages n L =-m=V/2 > 0. For |fi - V/2\ < V, 
and for V » D > |Q - V/2\,T, Eq. © gives, 



9d(D) « ~ 



A/2' 



i-e 



v 

TV 



^xA/2- 



The resummed T-matrix is then given by, 
r («) 



— 7TV Im 



1 R 
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(13) 
(14) 
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Expanding Eq.JTSlin powers of the interaction g, we have 
checked that this is consistent with Eqn. ©. 

For V < £1 < Tk the T-matrix is described by a cusp 
with power A. Apparently, for T < \Cl - V/2\ <C V this 
exponent changes to A/2 and one may want to interpret 
this as a new nonequilibrium scaling regime governed by a 
2K channel, instead of K channel behavior. However, de- 
coherence is so strong that this regime is 'unmeasurable': 
the extra factor in Eq. flU) 1 > {D/V) A / 2 > {T/V) A / 2 « 
e -( lnK )/ K always remains close to 1. Moreover, the peak 
at Q=V/2 is not much higher than the minimum at O=0, 

\m[T{V/2) - T(0)]/Im[T(0)] ~ (^) A ^ « 1. These 
results are shown schematically in Fig. [T] where we have 
used max{a, 6)=v a 2 + b 2 . 

Above, we have shown that a 2K channel fixed point 
cannot be stabilized by a finite voltage in an extended 
regime. However, if one considers a model where initially 
9lr -C gd such a regime becomes accessible. Remarkably, 
one can even calculate the none quili brium conductance 
exactly in this regime for K=l 16[. It has also been 



suggested by one of us [9] that a finite voltage stabilizes 
the 2K channel regime efficiently for K=l for sufficiently 
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small Qlr- It is therefore interesting to study this case 
also for K > 1. 

For gLR—0 the system is described by a 2K 
channel fixed point with a Kondo temperature 



T K =D {Kg d0 ) K e 



K p -l/g dQ 



and voltage has no effect as no 



currents can flow. However g^R is a relevant variable 
for V=0. While at the scale T K it is small, <7l_r(IV) ~ 
it grows below Tk with scaling dimension — 1/K 



9LRQ 

(compared to the exact result jj —\/{K + f)), thus in- 
ducing a flow back to the K channel hxed point which is 

. K 

reached at the scale T\ 
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New physics can arise in the regime T K < V < Tk 
governed by the 2i4T-channel fixed point. For exam- 
ple, the conductance in this regime is given by G m ~ 

j[gLR(V) ~ ( ) which crosses smoothly over 

to Eqn. (fTTj) for V ~ T£. Similarly, the decoherence 
rate in this regime is given by T ~ KVgLR(V) 2 ~ 

I J . To decide whether the range of valid- 
ity of the 2K channel regime is enhanced compared to 
the equilibrium case we compare T with T K and V for 
T K > V > Tt. We find 



T. 



K 



1 

K 



V_ 

T 



K 



1-2/ K 



r 

v 



i 

K 



V_ 

T~ 



K 



-2/K 



(16) 



While the ratio T/V is reduced in this regime, the de- 
coherence rate is always larger than T K for V/T K > K 
and approaches T/V s» 1/K as before for smaller volt- 
ages (up to small corrections of 0(\nK/K)). This implies 
that voltage bias does not enhance the regime where 2K 
channel physics is observable in contrast to the sugges- 
tion of Ref. || . 

Leads with different voltage differences: We now 
discuss the case where (for gLR=gd) a fraction p and 
hence Kp leads are at voltage V <C Tk, while K(l — p) 
are at V = 0. (In @| one has p = 1/2 and X = 2). The 
question arises whether the decoherence T again prohibits 
the flow to new fixed points for finite V. Analyzing as 
above the logarithmic corrections to the T-matrix in the 
regime D < V we obtain 



dgo 
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dlnD 


9a - 


dgdv 
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dlnD 
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go{(i-p)g 2 o+p(2g 2 dV )} 
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g d v{(l-p)g 2 +p(2g 2 dV )} 
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where gdv and go = 2gd are the coupling constants for the 
leads with and without an applied bias voltage respec- 
tively. The initial values gdv(V)=go{V)/2=g(V) /2 are 
obtained from Eq. ([9]) or (p~0|) . For our initial condition 
gdv < go, the solution apparently flows towards the fixed 
point g dv =Q,gQ=2/(K(l — p)) where only the channels 
without bias voltage contribute to screening, and pow- 
erlaws are governed by the exponent /S. p =2/{K(1 — p)). 



This is however misleading because the decoherence stops 
the flow towards this fixed point for all p. To see this note 
that close to the K channel fixed point gdv—go/^=^/K 
the RG equations are « ~2p/K 2 and « 

(1 —p)/K 2 . Therefore all changes arising from the flow 
from D = VtoD = T~ irpV/K remain for arbitrary 
p smaller than ^hi(V/T) < ^\n{K 2 ) < 1/K. There- 
fore the new fixed points can never be approached. 

In summary we have studied the nonequilibrium over- 
screened Kondo problem in the perturbatively solvable 
limit of a large number of leads K ^> 1. For this 
model, current induced decoherence is very small (T/V — 
Oil/K) -c 1). However, in this limit the renormaliza- 
tion group flows also become very slow. Our calculations 
show that the net result is that the flow to any new volt- 
age induced fixed points is stopped very effectively by 
the decoherence. We have made predictions for the split- 
ting of the T-matrix for V Tk which can be observed 
experimentally. An important open question is to study 
this nonequilibrium problem in the presence of an exter- 
nal magnetic field. 
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